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ABSTRACT 

We investigate the formation of protostellar clusters during the collapse of dense molec- 
ular cloud cores with a focus on the evolution of potential and kinetic energy, the 
degree of substructure, and the early phase of mass segregation. Our study is based 
on a series of hydrodynamic simulations of dense cores, where we vary the initial den- 
sity profile and the initial turbulent velocity. In the three-dimensional adaptive mesh 
refinement simulations, we follow the dynamical formation of filaments and protostars 
until a star formation efficiency of 20%. Despite the different initial configurations, 
the global ensemble of all protostars in a setup shows a similar energy evolution and 
forms sub-virial clusters with an energy ratio £'kin/|£'pot| ^ 0.2. Concentrating on 
the innermost central region, the clusters show a roughly virialised energy balance. 
However, the region of virial balance only covers the innermost ~ 10 — 30% of all the 
protostars. In all simulations with multiple protostars, the total kinetic energy of the 
protostars is higher than the kinetic energy of the gas cloud, although the protostars 
only contain 20% of the total mass. The clusters vary significantly in size, mass, and 
number of protostars, and show different degrees of substructure and mass segrega- 
tion. Flat density profiles and compressive turbulent modes produce more subclusters 
then centrally concentrated profiles and solenoidal turbulence. We find that dynami- 
cal relaxation and hence dynamical mass segregation is very efficient in all cases from 
the very beginning of the nascent cluster, i.e., during a phase when protostars are 
constantly forming and accreting. 

Key words: hydrodynamics - instabilities - stellar dynamics - stars: formation - 
stars: kinematics - turbulence 



1 INTRODUCTION 

In the current paradigm of star formation, most of the stars 

Con- 



form in a clustered environment (Lada et al 



cerning massive stars, studies by de Wit et al 



20031 



(2004 20051 



give an upper limit of 4 ± 2% of O-stars which could not 
be traced back to star clusters and which are thus candi- 
dates for the isolated formation of massive stars. Further 



work by Schilbach & Roser ( 2008 1 , Gvaramadze & Bomans 



(20081, and Pflamm-Altenburg fc Kroupa (20101 even allow 



for a lower fraction of O-stars that have to form in a clus- 
tered environment. Consequently, an understanding of star 
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formation is ultimately linked to the formation of clusters 
and stars within them. 

Over a huge spatial range of astrophysical objects and 
thus also during the collapse of a molecular cloud and the 
formation of a stellar cluster, the observed kinetic energy 
shows a robust scaling with the size of the object (|L ar-| 
son|1981| [Solomon et al.|1987|[Oss cnkopf & Mac Low 2002[ 
Heyer fc Brunt||2004| |Hily-Blant et al.„2008, ,Roman-Duvall 
et al.|20lT l. This global analysis, however, does not take into 
account the spatial and dynamical substructure of small- 
scale collapsing regions with sizes below 0.1 pc. The energy 
balance and virial state of the star-forming region may vary 
during the formation of the cluster and for different degrees 
of substructure in a cloud. Local changes in the dynamics 
may lead to different formation modes of the cluster and 
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alter the stellar distribution and the accretion process in a 
nascent cluster. 

Within a cluster, the distribution of stars is generally 
not uniform, but shows signatures of mass segregation with 
a tendency of more massive stars to be located closer to 
the centre of the cluster. This phenomenon is observed in 
many young clusters ([Hillenbrand fc Hartma nn|1998| 



Stoltel 



20M 



Espinoza ct al. 2009, Bontemps et al.^,2010, ^Gennaro et al.| 
2011). However, the detailed definitions of mass segregation 



and the regions where they apply lead to controversies. [Kirk] 
& Myers (20111 find mass segregation in small groups in 
Taurus, [Parker et al.| ( |20TT| ) find more massive stars to be 
inversely mass segregated, concentrating on the complex as 
a whole. In addition, there is a strong debate, whether the 
observed mass segregation in young clusters is primordial 
or due to dynamical interactions via two-body relaxation. 
One fundamental problem of that debate lies within the 
definition of mass segregation and the methods and tools 



to determine the segregation state. Allison et al. ( 2009 1 de- 



fine mass segregation as massive stars located close to other 
massive stars. |Kirk fc MyerS| ( ^2011J base their mass segrega- 
tion on the distance of the more massive stars to the centre 



of the local association. Maschberger & Clarke (20111 in- 



vestigated mass segregation of a collection of smaller cores 



and modified the model by Allison et al. ( 2009 1 to be more 



robust in case of outlier stars. In addition, they also used 
local surface density as a measure of mass segregation. Gen- 
erally, the substructure of the region in question plays a 
significant role in the explanation of the origin of mass seg- 
regation. Whereas the global system might not have enough 
time to dynamically relax, the small individual subclusters 
might well be able to reach a relaxed segregated state. In 
addition, the final mass segregation may crucially depend 
on how much degree of mass segregation is preserved during 
the merger of small subclusters, i.e. how much mass segrega- 
tion the merged structure can inherit from its constituents. 
Consequently, a combined investigation of the degree of sub- 
structure as a function of time, the energetic state of the 
cloud, the formation mode of stars within the cluster, and 
the formation of the clusters themselves is absolutely crucial 
to understand the mass segregation process. 

In this study we analyse the dynamical evolution of col- 
lapsing cloud cores and their virial state before and during 
the formation of protostars. In addition, we investigate the 
resulting substructure during the collapse and the possible 
degree of dynamical mass segregation for dense collapsing 
cloud cores in numerical simulations. We vary the initial den- 
sity profile as well as the initially imposed turbulent motions 
and analyse their impact on the later cluster structure. The 



simulations, which are taken from Girichidis et al. {2011al, 
hereafter Paper I, follow the collapse of the core and the 
formation of protostars. We find that the initial conditions 
have a large impact on the degree of substructure in a clus- 
ter and that the clusters show strong dynamical interactions 
between the protostars. As a result, the individual subclus- 
ters are very likely to have enough time for dynamical mass 
segregation. In contrast, for the global cloud, the time scales 
for dynamical relaxation are too long in comparison to the 
time scale at which stars form in these dense cores. Due 
to the strong dynamical interactions in the central region 
of the (sub)clusters from their formation onwards, it is ba- 



Table 1. Physical parameters of all setups 
Parameter Value 



cloud radius 


Ro 


3 X 10^^ cm 0.097pc 


total cloud mass 


Mtot 


100 Mq 


mean mass density 


(P> 


1.76 X 10"^* g cm-3 


mean number density 


(n> 


4.60 X 10^ cm-^* 


mean molecular weight 




2.3 


temperature 


T 


20 K 


sound speed 


Cs 


0.27 km s-l 


rms Mach number 


M 


3.28 - 3.64 


mean free-fall time 


tfS 


5.02 X 10* yr 


sound crossing time 


^SC 


7.10 X 10^ yr 


turbulent crossing time 


itc 


1.95 - 2.16 X 10^ yr 


Jeans length 


Aj 


9.26 X 10^ AU ft! 0.23 Ro 


Jeans volume 


V} 


1.39 X lO^i cm^ 


Jeans mass 


Mj 


1.23 Mq 



sically impossible to define primordial mass segregation in 
the simulated cores. 

The paper is organised as follows. Section [2| covers the 
description of the simulations and the key properties of the 
numerical setups. In section [3] we introduce the methods 
that we use to analyse the energy state, the degree of sub- 
structure, and the mass segregation. Section [4| presents our 
results, separately for the global cloud and the central or 
main subclusters. Section [S] and [6| comprise the discussion of 
the obtained results and the conclusions, respectively. 



2 NUMERICAL METHODS & INITIAL 
CONDITIONS 

The simulation data used in this study are the same as in 
Paper I, where a detailed description of the initial setups 
can be found. Here we only summarise the key parameters. 



2.1 Global Simulation Parameters 

We simulate the collapse of a spherically symmetric molec- 
ular cloud with a radius of i? = 0.1 pc and a total 
mass of 100 Mq. The resulting average density is (p) = 
1.76 X 10~^* g cm"'' and the corresponding free-fall time 
gives tff = 5.02 x 10* yr. The gas with a mean molecu- 
lar weight of /I = 2.3 is assumed to be isothermal at a 
temperature of 20 K, yielding a constant sound speed of 
Cs ~ 2.68 X 10* cm s~^. The Jeans length, Aj, and the corre- 
sponding Jeans mass Mj, calculated as a sphere with diam- 
eter Aj, are Aj = 9300 AU and Mj = 1.23 Mq, respectively. 
Table [l] provides an overview of all physical parameters. 



2.2 Numerical Code 

The simulations were carried out with the astrophysical 
code FLASH Version 2.5 ( [FryxeU et aL][2000[ ). To integrate 
the hydrodynamic equations, we use the picccwise-parabolic 



method (PPM) by iColella fc Woodward ( 1984 ) . The compu- 
tational domain is subdivided into blocks containing a fixed 
number of cells with an adaptive mesh refinement (AMR) 
technique based on the PARAMESH library ( [Olson et al.[ 
1999|). 
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Table 2. Numerical simulation parameters 



Parameter 




Value 


simulation box size 


^box 


0.26 pc 


smallest cell size 


Ax 


13.06 AU 


Jeans length resolution 




^ 8 (6*) cells 


max. gas density 


Pmax 


2.46 X 10-1"' g cm-3 


max. number density 




6.45 X 10^ cm-3 


sink particle accretion radius 




39.17 AU 



* at highest level of refinement 

2.3 Resolution and Sink Particles 

The simulations were run with a maximum effective resolu- 
tion of 4096'' grid cells, corresponding to a smallest cell size 
of Ax « 13 AU. In order to avoid artificial fragmentation, 
the Jeans length has to be resolved with at least 4 grid cells 
(Truelove et al. 19971. To resolve turbulence on the Jeans 



scale, however, a significantly higher number of cells is re- 
quired. Fcdcrrath et al. (|2011|) find a minimum resolution of 
about 30 cells per Jeans length. Due to the high computa- 
tional demand, we only use 8 cells in the current runs, so we 
likely miss some turbulent energy in our cores, which pro- 
vides additional support against gravitational collapse. We 
might thus slightly overestimate the amount of fragmenta- 
tion and imderestimate the formation times of protostars. 
It must be noted, however, that this is a general limitation 
of all present star cluster formation calculations because re- 
solving the Jeans length with more than 10 — 20 cells can be 
computationally prohibitive. Additionally, in order to termi- 
nate local runaway collapse in a controlled way, we use sink 



particles (see e.g., Bate et al. 1995 Krumholz et al. 2004 



Federrath et al.|2010a |. They are introduced at the highest 



level of the AMR hierarchy. A necessary but not sufficient 
criterion for the formation of sink particles is that the gas 
density needs to be higher than the threshold value 



Pmax 



AG {3 Ax) 



2.46 X 10""g cm"^ 



(1) 



If a cell exceeds this density, a spherical control volume with 
a radius of 3Aa; is investigated for the following gravitational 
collapse indicators ( Federrath et al.|2010a k 



The gas is converging along all principal axis, x, y, and 



z, 



• has a central minimum of the gravitational potential, 

• is Jeans-unstable, 

• is gravitationally bound, and 

• is not within the accretion radius of an already existing 
sink particle. 

If the collapse criteria are fulfilled, an accreting Lagrangian 
sink particle is formed. This sink particle is then identified 



1995 Wuchterl & 



as an individual protostar ( Bate et al 

Klessen 20011. Table |2] lists the simulation and resolution 
parameters. 



(iii) Power-law profile, p oc r 

(iv) Power-law profile, p cx r 



-1-5 (PL15) 
° (PL20). 



A detailed description of the profiles can be found in Paper I. 

The turbulence is modelled with an initial random ve- 
locity field, originally created in Fourier space, and trans- 
formed back into real space. The power spectrum of the 
modes is given by a power-law function in wavenumber space 
(k space) with _Ek oc fc~^, corresponding to Burgers turbu- 
lence, consistent with the observed spectrum of interstellar 



turbulence (e.g., Larson|l981 Ossenkopf fc Mac Low|2002 



Heyer & Brunt 20041. The velocity field is dominated by 
large-scale modes due to the steep power-law exponent, —2, 
with the largest mode corresponding to the size of the sim- 
ulation box. Concerning the nature of the k modes, com- 
pressive (curl-free) modes are distinguished from solenoidal 
(divergence-free) ones. The simulation uses three types of 
initial fields: purely compressive fields (c), purely solenoidal 
(s), and a natural (random) mixture (m) of both. The choice 
of these different turbulent fields was motivated by the 
strong impact of the nature of the modes on the cloud evolu- 



tion, found by Federrath et al. (2008 2010b I. Note however 



that only decaying turbulence with compressive, solenoidal, 
and mixed modes are considered here. 

All setups have supersonic velocities with an rms Mach 
number M = «rms/cs ranging from M = 3.28 — 3.54 with 
an average of (M) — 3.44. The sound crossing time through 
the entire cloud is tsc{Ro) ~ 7.10 x 10^ yr, and the time for 
gas with an average velocity of (A^) Cs to cross the cloud is 
itc(-Ro) = 2.06 x 10^ yr, respectively. 

We combine four density profiles with six different tur- 
bulent velocity fields (three different compositions of modes 
with two different random seeds each). Table [s] shows a list 
of all models. 



3 CLOUD AND CLUSTER ANALYSIS 

In this section we briefiy motivate and summarise the meth- 
ods we used to analyse our simulation data. 



3.1 Energy Analysis 

The global energy partitioning of a gas cloud can be 
quantified by the ratio of kinetic to the potential energy 
Ei^in/\Epot\, where a value of 0.5 corresponds to a virialised 
cloud. During the collapse of the cloud and the collapse of 
fragments into protostars, potential energy is converted into 
kinetic energy and transfered from the smooth gas to rela- 
tively compact protostars. In order to investigate the energy 
evolution of the collapse, we analyse the energy budget for 
the gas and the protostars separately. 

The total kinetic energy of the gas is calculated by sim- 
ply summing over all cells in the cloud 



P — 1 / 2 I 2 I 2 

-L/kin.gas — ^ / ^ y'i^i^x ' ^i^y ^ Vi^^ 



(2) 



2.4 Initial Conditions 

The following four density profiles were used: 

(i) Top-hat profile, p = const (TH) 

(ii) Rcscaled Bonnor-Ebert sphere. (BE) 



The kinetic energy of the protostars, -Ekin,sink, is found anal- 
ogously. For the potential energy of the gas we integrate nu- 
merically over radial bins around the centre of mass, yielding 



^pot.gas 



(r) 



GM{r)Am{r) 



(3) 
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Table 3. List of the runs and their main properties 



name M total total tsim tsim/tff ^sink Mnax nf°^ (s) {s)^ {m)^ Q 



TH-m-1 


3.3 


0.075 


0.047 


48.0 


0.96 


311 


0.86 


5.50 




10* 


6.51 


0.42 


0.11 


0.26 


TH-m-2 


3.6 


0.090 


0.047 


45.5 


0.91 


429 


0.74 


8.00 


X 


10* 


8.51 


0.65 


0.14 


0.21 


BE-c-1 


3.3 


0.058 


0.039 


27.5 


0.55 


305 


0.94 


1.70 


X 


10^ 


3.11 


0.16 


0.09 


0.53 


BE-c-2 


3.6 


0.073 


0.039 


27.5 


0.55 


331 


0.97 


3.60 


X 


10" 


5.68 


0.31 


0.08 


0.27 


BE-m-1 


3.3 


0.053 


0.039 


30.1 


0.60 


195 


1.42 


3.20 


X 


10^ 


1.10 


0.13 


0.13 


1.03 


BE-m-2 


3.6 


0.074 


0.039 


31.9 


0.64 


302 


0.54 


2.48 


X 


10*^ 


1.46 


0.13 


0.09 


0.74 


BE-s-1 


3.3 


0.055 


0.039 


30.9 


0.62 


234 


1.14 


3.70 


X 


10^ 


0.52 


0.11 


0.14 


1.30 


BE-s-2 


3.5 


0.074 


0.039 


35.9 


0.72 


325 


0.51 


3.20 


X 


10^ 


1.43 


0.21 


0.14 


0.68 


PL15-C-1 


3.3 


0.056 


0.038 


25.7 


0.51 


194 


8.89 


2.42 


X 


10*^ 


1.99 


0.11 


0.06 


0.71 


PL15-C-2 


3.6 


0.068 


0.038 


25.8 


0.52 


161 


12.3 


1.66 


X 


10* 


7.82 


0.45 


0.09 


0.21 


PL15-m-l 


3.3 


0.050 


0.038 


23.8 


0.48 


1 


20.0 
















PL15-m-2 


3.6 


0.071 


0.038 


31.1 


0.62 


308 


6.88 


2.66 


X 


10*^ 


1.21 


0.11 


0.11 


0.99 


PL15-S-1 


3.3 


0.053 


0.038 


24.9 


0.50 


1 


20.0 
















PL15-S-2 


3.5 


0.069 


0.038 


36.0 


0.72 


422 


4.50 


1.11 


X 


10^ 


1.01 


0.16 


0.19 


1.20 


PL20-C-1 


3.3 


0.042 


0.029 


10.7 


0.21 


1 


20.0 

















The acronym for the run is shown in the first column, where the first part indicates the density profile, the middle letter the turbulent 
mode ('c' for compressive modes, 's' for solenoidal modes and 'm' for a natural mix of both), and the number at the end of each name 
the random seed for the turbulence. The initial energetic state is given by the Mach number M, and the ratios of kinetic and thermal 
energy to the potential energy, t^i^ and tsim/*ft show the simulation time, N^i^]^ the total number of protostars, and Mmax the mass of 
the most massive protostar. The stellar number density is shown in column nf^°^ . The global cluster properties are given as the mean 
separation between the protostars (s), the normalised mean separation (s)^, the normalised mean length of the minimal spanning tree 

(m)^, and the ratio Q. 



where G is Newton's constant, A/(r) the enclosed mass in- 
side radius r, and dm(r) the mass in the radial shell with 
thickness dr. The potential energy of the protostars can be 
calculated by summing over the point masses 

rrii rrij 
■ rA 



(4) 



However, in order to avoid the formation of hard binary 
systems and resulting very small time steps, we apply a 
softening term in the computation of the gravitational force 
between the protostars. For the softening we use the energy- 



conserving formalism described in Price & Monaghan ( 2007 1 
which yields a potential energy of 



E, 



pot, sink 



(5) 



with a kernel function ip (see Appendix|A]| . On the one hand, 
the applied softening artificially prevents the formation of 
hard binaries and close orbits of particles in the simulation. 
On the other hand, it is questionable to what extend hard 
binaries can form in the early evolutionary phase. Numeri- 
cally, the protostars are point objects with arbitrarily close 
separations. Physically, the protostars are very young and 
still in the contraction phase. Consequently, they have a rel- 
atively large sizes and low density contrasts in comparison 
to main sequence stars. Therefore, a dynamical treatment as 
extended gas spheres might well be more realistic. However, 
the detailed substructure inside the sink particle radius and 
the resulting dynamics is not captured in our simulations. 

The internal motions of the gas and the protostars are 
quantified using the mass-weighted velocity dispersion 



(6) 



where k € {x, y, z} and (uk) is the mean velocity in dimen- 



sion fc, 



(Uk) = 



J2i '^k,: 



(7) 



The three-dimensional velocity dispersion is then given by 



0"3D 



(8) 



In the simulations we calculate ctsd using each component 
of the velocity. For the one-dimensional velocity dispersion 
we assume the same value for all three components and thus 
use (TiD ~ f7"3D/\/3. So far we have only considered the tur- 
bulent contribution to the velocity dispersion. Including the 
thermal contribution, the total dispersion along the line of 
sight is 



cy tot 



(9) 



3.2 Subclustering 



Depending on the interplay between turbulent motions and 
the central collapse of a cloud, the spatial distribution of 
protostars may vary significantly (see Paper I). In order to 
analyse the clustering properties of our protostars, we use 
the Q value ( |Cartwright fc Whitworth|2004 1 



Q 



(10) 



of the clusters. Here, {s)^ is the normalised mean separation 
of the protostars and {tti)^ is the normalised mean length 
of the edges of the minimal spanning tree (MST), where 
the edge is the distance between two protostars. For a de- 
tailed discussion of the motivation for this definition of Q 



see Cartwright & Whitworth ( 2004 1 
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The distribution function p(s) describes the probability 
of two protostars to be separated by the distance, s. We 
discretise p{s) with A'^bin bins for the entire cluster, leading 
to an equal-sized bin width of As = 2_Rc/-^bin, where Rc is 
the cluster radius. The normalised number of pairs in bin i 
can thus be expressed with 

^«=iVc(iVc-l)A.- (11) 

Here Ni denotes the number of pairs with a distance in the 
range [iAs, {i + l)As) and Nc{Nc — l)/2 is the total num- 
ber of separations for A^c cluster members. Multiple peaks in 
the distribution function are related to subcluster structure, 
which gives higher counts at low distances due to the small 
separations within each subcluster and higher counts at a 
larger separation due to the large distance between the sub- 
clusters. In case of no distance degeneracy between subclus- 
ters, the number of peaks equals the number of subclusters. 
The mean value (s) of all Nc{Nc — 1) particle separations 

gives a measure for the mean distance between particles in 
the set. 



The MST is calculated usin g thelGower fc Ross| ( |1969 l 
description of Prim's algorithm (Prim|1957l. The more par- 



ticles are confined in an observed area, the smaller is the 
mean edge of the tree. The resulting decrease of the mean 
edge due to the increasing number of nodes in the tree has 
to be corrected by a dimensionality factor. The correction 
factor for the three-dimensional cluster model with cluster 
volume V was set to 



Nc-l 



(13) 



taken from |Schmeja fc Klessen] ( |2006[ ). 

For stellar clusters with a smooth radial density gra- 
dient, Q ranges from 0.8 — 1.5, corresponding to a radial 
density distribution of particles n oc with = to 2.9. 
Clusters with substructure have Q = 0.8 — 0.45, decreasing 
with increasing degree of subclustering. A detailed relation 
between Q, rj and the degree of subclustering can be found 



in Cartwright & Whitworth (20041. 



3.3 Mass Segregation 

A set of stars or protostellar objects may show a mass- 
dependent spatial distribution within a cluster. In a mass- 
segregated cluster, massive objects tend to be located closer 
to the centre of the cluster, whilst low-mass objects occupy 
regions of larger radii. We quantify the degree of mass seg- 



regation using the MST as described in | Allison et al. 
with the mass segregation ratio (MSR) 



(20091 



A^ 



(^iiorm) _|_ *^norm 



massive 



(14) 



massive 



The ratio describes, how large the spatial spread of the most 
massive stars is, compared to the spatial spread of a random 
choice of stars. How many most massive stars are counted 
and compared to an equal amount of random stars should 
not be fixed, but rather treated as a free parameter, which 



we name A''mst. In order for the MST of the random set of 
stars to be a good measure for the average spread, we need to 
pick many sets of random stars and average over the individ- 
ual lengths of the MST. We set the number of sets to 500 as 



suggested by Allison et al. (20091. With the average length 
(Znorm) of thcse 500 scts and the length of the A^mst most 
massive stars, /massive, we then determine the degree of mass 
segregation. The error is computed with the standard devia- 
tion (Tnorm of (/norm). If Amsr takes values siguificautly larger 
than unity, the A^mst most massive stars are located much 
closer to one another than the same amount of randomly 
picked stars. Hence the system shows mass segregation. In 
the opposite case (Amsr ^ 1) the most massive stars have 
much larger distances between one another than a set of 
random stars in the cluster and the system shows inverse 
mass segregation. A^mst is basically a free parameter that 
we loop over starting from 2 up to half of the total number 
of sink particles, in order to determine the number A''mst up 
to which the system is mass segregated, i.e., Amsr > 1. The 
weak point of this method is its sensitivity to massive outlier 
protostars, in particular protostars among the A^^mst most 
massive stars that are clearly located between subclusters in 



an environment of many separate clusters (see, Maschberger 



& Clarke 2011 1. As we do not investigate the total cloud with 



this method but only reduced clusters without outliers (see 
section 



4.4 1 our results are not affected by this behaviour. 



Mass segregation can either originate from dynamical 
Ai'-body relaxation or is primordial in nature, where the lat- 
ter case means the more massive stars form closer to the 
centre. In order to analyse whether mass segregation is pri- 
mordial or due to dynamical processes, we use the mass 
segregation time ( |Spitzer^l969j ), 



tseg(M) 



(m) 
M 



(15) 



with (m) being the average mass of all stars in the cluster 
and M the mass of the star in question. The relaxation time 
iroiax can be expressed in terms of the number of stars TV, the 
radius of the cluster Rc, and the stellar velocity dispersion 
a, yielding for the mass segregation time (e.g. [Binney "fc| 
Tremaine||1987[ ), 

m) TV Rc 
M 



tscg(TW") 



(16) 



81nTV a ' 

By setting the time according to different stages in the sim- 
ulation, one can obtain the minimum mass down to which 
stars had enough time to dynamically mass segregate. Care 
must be taken when applying the mass segregation time to 
hydrodynamic collapse simulations. In contrast to old stel- 
lar clusters, where there is no or very little interstellar gas 
left and consequently TV, (m) and M do not vary with time, 
hydrodynamic collapse simulations follow the formation of 
protostars from the beginning of the collapse. Not only do 
protostars form at different times, they also accrete further 
gas from the surrounding dense medium in which they were 
born and are subject to gas drag forces. The number of pro- 
tostellar objects TV, their individual masses M, their mean 
mass (m), and the cluster radius R are therefore strongly 
varying with time. Consequently, the mass segregation and 
the minimum segregated mass for a given time can not be 
calculated for the total set of objects as a whole. Instead, the 
possibility of being segregated within the cluster has to be 
estimated for each star individually by taking into account 
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the formation time and the growing mass of the star due to 
accretion. 



4 RESULTS 
4.1 Overview 

We follow the cloud collapse until 20% of the mass is ac- 
creted by sink particles. The simulation time, the number of 
formed protostars, the mass of the most massive protostar 
and the key parameters of the substructure of the cluster 
are listed in table [3] A column density plot at the end of 
each simulation is shown in Paper I, figures 4 and 5. 

The TH profile takes the longest time to form gravi- 
tationally collapsing regions and to capture 20 Mq in sink 
particles. During this time, approximately 45 — 50 kyr, the 
turbulent motions can compress the gas in locally discon- 
nected areas, leading to distinct subclusters of sink parti- 
cles. The stronger mass concentration in the centre of the 
BE setups and the resulting shorter collapse and sink parti- 
cle formation time suppresses the formation of disconnected 
subclusters in favour of one main central cluster (see mor- 
phology in Paper I). The corresponding PL 15 profiles show 
a very similar overall cloud structure to the BE runs, but 
significantly different stellar properties. Due to the much 
stronger gas concentration in the centre of the cloud, all 
PL15 setups form a protostar very early in the simulation. 
This initial central protostar accretes the surrounding gas at 
a high rate and can grow to a massive protostar before the 
turbulent motions eventually form collapsing filaments and 
trigger fragmentation. The PL15 setups with turbulent fields 
m-1 and s-1 (PL15-m-l, PL15-S-1) do form dense filaments, 
but no further sink particles until the first protostar reaches 
a mass of 20 Mq. In case of multiple sink particles, the clus- 
ters are more compact than in the corresponding BE case. 
The PL20 profile only forms one single sink particle due to 
the very strong mass concentration. The central protostar 
forms very early and accretes gas at an almost constant rate 
of ~ 2 X 10~^ Mq yr~^, close to the analytical value of a 



10 



highly unstable singular isothermal sphere (Shu 1977 Pa- 
per I). This results in a total simulation time of only 11 kyr, 
which is not enough for turbulent motions to form filaments 
and further sink particles. 

The following discussion of the cluster properties there- 
fore abstains from a detailed description of the setups PL15- 
m-1, PL15-S-1 and PL20-C-1. 



4.2 Energy Evolution of the Global Cloud 

In order to better understand the energy evolution, we sepa- 
rately analyse the gas and sink particle contributions to the 
total energy. 

All setups are gravitationally very unstable and start 
to collapse immediately. As a result, the initial random ve- 
locities of the gas are reoriented towards the direction of 
the central acceleration. The total kinetic energy strongly 
increases with time due to the infall motion. Figure [l] shows 
a representative example of the kinetic over the potential 
energy of the gas as a function of radius for different times 
in the simulation. The cloud starts in a strongly sub-virial 
state and exceeds a ratio of kinetic to gravitational energy of 
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Figure 1. Ratio of kinetic over potential energy of the gas for 
the BE-m-1 setup as a function of radius for different times. The 
cloud evolves from a strongly gravitationally dominated state to 
an energy state with i?kin,gas/|£'pot,gas| > 0.5 at the end of the 
simulation. 



0.5 for the entire cloud after roughly 20 kyr. Within a radius 
of 10^ AU the ratio reaches values greater than unity and 
diverges in the very central region. This behaviour can be 
explained by a simple estimate using a singular isothermal 
sphere, which is characterised by an initial density profile 
p oc and approaches a free-fall density profile p oc r~^^^ 



Whit- 



field 



inside the head of the rarefaction wave (Shu 1977 
I worth fc Summers||1985[ ). The corresponding velocity 
scales as w oc r~^' . The resulting potential energy scales 
as -Epot oc r^, while the kinetic energy follows a relation 
-Ekin oc r^/^. Consequently, the ratio -Ekin/|i?pot| scales as 
r^^^'^ and diverges for small radii, indicating that the inner- 
most part of the cloud is dominated by kinetic energy. 

The different initial density profiles as well as the dif- 
ferent formation modes of protostars lead to different ra- 
dial distributions during the collapse. A comparison of 
i?kin,gas/|£'pot,gas| for all sctups at the end of the simulation 
is shown in figure|2] A significant difference is found between 
the simulations with only one protostar (dotted lines) and 
the ones that form many protostars (solid lines). The three 
setups with only one protostar show much higher values for 
most of the cloud and a steeper slope. This is not surpris- 
ing because the gas in the central region can fall towards 
the central particle without being disturbed by other sink 
particles and their A^'-body interactions. In case of multiple 



protostars the ratio Skin 



pot, gas 



shows a large scatter 



close to the central region (-R < 4 x 10'^ AU), which can 
be explained by the local variations in the sink particle po- 
sitions and motions, and the resulting impact on the gas. 
The scatter in the energy ratio is significantly lower in the 
outskirts of the cluster. 

The average value as well as the spread of 
i?km,tot/|i5pot,tot| increase when the sink particles' mass is 
included in the virial analysis (see figure |3|. There is no sys- 
tematic correlation between the various initial conditions 
and the ratio of the energies. The fact that including the 
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Figure 2. Ratio of kinetic over potential energy of the gas for 
all profiles at the end of the simulation at 20% star formation 
efficiency. The dotted lines indicate the runs with only one proto- 
star, the runs with multiple protostars are shown with solid lines. 
Note that the physical times differ strongly between 11 and 48 
kyr for the different setups, see table [3] 



protostars leads to higher values, shows that the cluster con- 
tributes more to the kinetic rather than the potential energy. 
At this point, we want to emphasise that the computation of 
the potential energy with and without gravitational soften- 
ing shows different values that vary by a factor of a few. This 
also influences the kinetic energy evolution of protostars in 
close encounters. Considering the large values of kinetic over 
potential energy and the large spread in the central regions 
of the clouds, this does not affect the overall result that 
the central region is strongly dominated by kinetic energy. 
On longer evolutionary time scales, one has to take into ac- 
count that the very young protostars are still large objects 
that slowly contract. The effects of protostellar collisions at 
different stages, the resulting mergers and formed binaries 
may have different effects on the global energy evolution. 

A comparison of the kinetic energy of the sink particles 
and the gas (-Ekin,sink/-Bkin,gas) is plotted in figure [i] The 
ratio is above unity for all simulations with many proto- 
stars (solid lines). Although the protostars account for only 
20% of the total mass at the end of the simulation, their 
kinetic energy dominates the total kinetic energy budget of 
the cloud. Again, the setups with only one protostar consti- 
tute an exception (dotted lines). In these cases, the kinetic 
energy of the protostar is significantly lower, which can be 
explained by accretion fiows from opposite directions that 
result in an almost vanishing net momentum transfer onto 
the protostar (see figure [4|. The dashed-dotted line shows 
£'kin,sink/£'kin,gas for the TH-m-2 setup. As the cloud in this 
run forms two distinct subclusters with a central void be- 
tween them (see right part of figure [9]l the total kinetic en- 
ergy of the few protostars between the subclusters is rela- 
tively low. 

As a link to observable properties of star-forming re- 
gions we calculate the velocity dispersion for the entire cloud 
as a function of time. Here we assume isotropy of the motions 
of the gas and restrict our analysis to the one-dimensional 
velocity dispersion ctid- Because of the initial random turbu- 
lence, the velocity dispersion of the gas shows anisotropies. 
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Figure 3. Same as figure [5] but for the total energies (protostars 
and gas). 
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Figure 4. Ratio of kinetic energy of the sink particles to the ki- 
netic energy of the gas for all profiles at the end of the simulation. 
The inner region {Rc < 4 AU) varies extremely because of the 
slight offsets of the centre of mass and the centre of the cluster 
and is not shown. 



which tend to reduce during the simulation. Initially, the 
deviation from isotropy is of the order of 10 — 20%. During 
the simulation the value shows variations but decreases to 
about half of the initial value {Aa/a ~ 5 — 10%), averaged 
over all simulations. There is no clear trend with the var- 
ied initial conditions and the number of protostars. FigurejS] 
shows the turbulent velocity dispersion ctid for the gas for 
all runs. Initially, a just refiects the initial turbulent velocity, 
the increasing values correspond to the additional infall mo- 
tion. The significantly lower values for the TH profiles are 
simply due to the delayed dominant central collapse. The 
formation of disconnected subclusters reduces the global in- 
fall speed in comparison to the other setups with one central 
cluster. The combined velocity dispersion for gas and sink 
particles can be seen for the TH profiles in figure [6] The 
plots for the other setups look similar. As shown in figure |4j 
the protostars contain a significant fraction of the kinetic 
energy. Therefore, the total value including sink particles is 
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Figure 5. Velocity dispersion for the gas for all runs. The values 
increase over time due to the increasing infall motion. 
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Figure 7. Ratio of kinetic to potential energy as a function of 
time for the sink particles only. The time was adjusted to the 
formation of the first sink particle in the setup. 
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Figure 6. Velocity dispersion for the TH setups. The lower curves 
correspond to the gas only, the higher curves include the turbulent 
contribution of the sink particles. As the sink particles contain a 
significant fraction of the kinetic energy, the curves including sink 
particles are remarkably higher. 



remarkably higher. None of the curves saturates during the 
simulated time, which can be explained by a simple free-fall 
approximation. The maximum speed that can be reached by 
free-falling gas is of the order of Ro/tg ~ 2 km s~^, where 
Ro is the cloud radius and the global free-fall time. None 
of the setups needs more than a free-fall time to convert 
20% of the gas mass into stars when we stop the simulation, 
so no setup had enough time to reach the limiting free-fall 
velocity dispersion of 2 km s~^. 

With a focus on the nascent cluster as an A''-body sys- 
tem, we also analyse the virial state of the sink particles 
without including the contributions of the surrounding gas. 
To do so, we treat the protostars as point masses and we 
calculate the gravitational potential via direct summation 
(equation [5| . The corresponding ratio of kinetic to poten- 
tial energy for the sink particles is shown in figure [T] ex- 
cluding the runs with only one protostar. The time axis in 
the plot is adjusted to the time when the first condensation 
was created. In the case of all PL15 profiles with multiple 
sink particles, the second and further sink particles formed 
with a large delay after the first sink particle. Therefore, the 
curves for the PL15 profiles start at times t — to > 10 kyr 



(see Ati2 in table [4|. The first protostars form with the ve- 
locity that the collapsing condensation inherits from the gas 
motion. The positions at which they form is determined by 
the structure of the filaments or the fragmenting disc. As 
they form independent from one another, sometimes even in 
separate subclusters, their kinetic and potential energies are 
uncorrelated. As the protostars form with the velocity of the 
parental gas cloud and because they are usually separated by 
a large distance, the initial values of i5kiii,sink/|i5pot,siiik| are 
very high. Soon after their formation, the protostars dynam- 
ically decouple from the gas and move towards the central 
region of the nascent cluster. The system begins to virialise, 
leading to decreasing values of i?kin,sink/|£'pot,sink|- Without 
the formation of subsequent protostars, the system would 
quickly reach a virialised state. However, as this process con- 
tinues, the energy ratio of the total cluster is influenced by 
the virial state of the newly formed objects. If they form at 
time ti at position with velocities Vi smaller than the virial 
velocity Uviriai(?'i, they lead to a decreasing energy ratio. 
A quick analytical estimate illustrates, why this behaviour 
is expected. The virial velocity is given by 



_ fGMci\ 



1/2 



(17) 



with the mass and radius of the cluster Mci and Rci- As a 
lower limit, we can assume a constant stellar density in the 
cluster over time, p*, which relates the cluster radius to the 
cluster mass like Rci{t) — (3Mci(t)/(47rp»))^^'' and thus the 

1/3 

virial velocity in this lower limit follows Wviriai.iow oc M|-,j , 
increasing with time as the total mass of the cluster in- 
creases. Of course, the velocity of the gas is also increasing 
over time due to the collapse of the cloud. However, as shown 
in figure [5] the velocity dispersion of the gas increases over 
time by a factor of only 3 at most. In addition, figure [6] 
illustrates that the kinetic contribution of the protostars is 
remarkably larger than that of the gas. In order for the lower 
limit virial velocity, Wviriai.iow, to be higher than the average 
gas velocity, the cluster mass must grow by a factor of 27 
during the entire simulation, which can be achieved. Con- 
sidering the fact, that the stellar density also increases, the 
virial velocity will be even higher. Consequently, the newly 
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Figure 8. Energy balance Skin^sink/|£^pot,siiik 
normalised number of protostars, Ni^i^i^/N^i^]^ tot - The protostars 
were sorted by their distance from the centre of the cluster. In all 
cases, only the innermost ~ 10 — 30% of the stars form a cluster 
with virial or super-virial energy balance. 



formed stars, which inherit the low gas velocity, tend to de- 
crease the energy ratio. 

The larger the cluster, the lower is the available mass 



in the central region of the cluster (see also, Girichidis et al. 
2011b I. Therefore, new protostars must form at increasingly 



larger radii. In order to show that these new stars are the 
ones that push the ratio -Ekin,sink/|i?pot,sink| to lower than 
virialised values, we calculate the ratio as a function of the 
fraction of total protostars. Figure [8] shows the energy ratio 
with the protostars sorted by their distance from the centre 
of the cluster. In all cases, only the innermost ~ 10 — 30% 
have a virial or super-virial energy balance. The majority of 
the nascent cluster has an overall sub- virial energy partition. 
But we expect that the ensemble virialises on a dynamical 
time scale as soon as star formation stops in the cluster 
region. 

From our simulations we conclude that a detailed en- 
ergy analysis can only be performed properly, if both pro- 
tostars and gas are included in a self-consistent way. In 



Table 5. Subcluster properties from the TH setups 

subclustcr AT^ink (s) [10^ AU] (m> [10^ AU] 



SCI 


67 


1.13 


0.31 


0.69 


SC2 


182 


0.49 


0.19 


1.36 


SC3 


232 


0.51 


0.16 


1.19 


SC4 


100 


0.23 


0.10 


1.27 



For each subcluster the number of sink particles N^i^]^, the mean 
separation (s), the mean MST length (m) and the Q value are 
shown. SCI shows signs of sub-structure indicated by a Q value 
slightly below the critical transition value of 0.8. SC2, SC3 and 
SC4 have values of Q > 1.2 which indicates a smooth internal 
structure. 



turn, the remaining gas is essential to the virial state of 
the nascent cluster. 



4.3 Global Cluster Properties 

In this section we discuss the spatial distribution of the pro- 
tostars in the simulated cloud. We begin with an analysis 
of all protostars in the simulations in order to measure the 
cluster properties of the cloud as a whole. A detailed inves- 
tigation of individual subclusters without outlier protostars 
is presented in section 1X4] and below. 



4.3.1 TH runs 

Both setups with initially uniform density distribution show 
distinct subclusters as illustrated in figure [9] We selected 
the four biggest subclusters for further analysis and named 
them SC1-SC4. The other subclusters have too few proto- 
stars for a statistical analysis. Note that subcluster SCI is 
not very compact in the centre. Therefore, our reduction 
algorithm does not exclude the outliers, which yields the 
relatively large radius. 

The distribution function of the separations between 



the particles as well as the Q-value (see equation 10 1 of the 
entire cloud is shown in figure [To] TH-m-1 shows three dif- 



ferent peaks in the distribution function (see equation 111: 
the one at 9,000 AU corresponds to the distance of SC2 
to SCI, the peak at 13,000 AU to the degenerate distance 
of SC2 to SC5 and SC6, and the last peak describes the 
distance from the upper subcluster SCI to SC5 and SC6, 
which is also degenerate within the width of the distance 
bin. TH-m-2 shows two main subclusters corresponding to 
the peak at 15,000 AU in the plot. The degenerate distance 
between SC3 and SC7 as well as SC4 and SC7 can be seen 
as small peak in the distribution at 13,000 AU. The Q value 
of the entire cloud shows strong variations at the beginning 
of stellar formation due to the different regions of the cloud 
where the sink particles are created. Having established the 
subclusters, Q shows roughly constant behaviour at a value 
of Q ~ 0.2 for both runs. 

The key properties for the subclusters SC1-SC4 are 
listed in table [5] The protostars in SCI have significantly 
larger mean separations between one another and a Q value 
of ~ 0.7, slightly lower than the threshold value to sub- 
structure of 0.8. The other three subclusters have very sim- 
ilar Q, indicating a smooth stellar distribution. 
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Table 4. Reduced cluster properties for the simulations with many sink particles 



Cluster A''sink Mq (M) _Rc K'"^ *1 ^*12 Cid t^^lax ^avail ^scg /sog SCg? 

[Mq] [Mq] [kAU] [pc-3] [kyr] [kyr] [km/s] [kyr] [t,,i,J 



SCI (TH-m 


1) 


67 


4 


2 





063 


2 


74 


6 


82 


X 


lO*' 


32 








008 


1 


86 


8 


01 


2 


00 


19 





28 





SC2 (TH-m 


1) 


182 


10 


4 





057 





97 


4 


18 


X 


10* 


29 


9 





776 


2 


90 


4 


02 


4 


31 


72 





40 




SC3 (TH-m 


2) 


232 


9 


4 





041 


1 


00 


4 


86 


X 


108 


26 


5 





709 


2 


82 


5 


17 


3 


53 


82 





35 


+ 


SC4 (TH-m 


2) 


100 


5 


7 





057 





45 


2 


30 


X 


10" 


28 


5 





933 


3 


00 


1 


12 


14 


31 


77 





77 


++ 


BE-c-1 




192 


11 


4 





060 





64 


1 


53 


X 


10" 


14 


9 





279 


3 


61 


2 


21 


5 


58 


81 





42 





BE-c-2 




275 


15 








055 


5 


05 


4 


47 


X 


lO'' 


15 


1 





764 


2 


43 


34 


93 





33 


8 





03 





BE-m-1 




99 


11 


9 





121 





24 


1 


50 


X 


lOio 


19 


6 





052 


5 


91 





30 


34 


66 


84 





85 




BE-m-2 




255 


15 


7 





061 


1 


39 


1 


99 


X 


10* 


20 


2 





086 


3 


20 


6 


81 


1 


71 


67 





26 


-1- 


BE-s-1 




190 


16 


1 





085 





50 


3 


18 


X 


10^ 


21 


5 





083 


4 


64 


1 


35 


6 


92 


100 





53 


++ 


BE-s-2 




288 


16 


7 





058 


2 


12 


6 


33 


X 


10^ 


22 


3 





004 


2 


85 


12 


97 


1 


05 


43 





15 




PL15-C-1 




170 


17 








100 


1 


46 


1 


14 


X 


10* 


1 


1 


13 


5 


4 


02 


4 


11 


2 


69 


37 





22 




PL15-C-2 




79 


14 


8 





187 


1 


64 


3 


75 


X 


10^ 


1 





15 


5 


2 


50 


1 


45 


6 


43 


13 





16 





PL15-m-2 




240 


15 


6 





065 


1 


00 


5 


03 


X 


10* 


1 





13 


3 


4 


47 


3 


34 


5 


03 


68 





28 




PL15-S-2 




396 


18 


5 





047 


1 


46 


2 


67 


X 


108 





9 


10 


3 


3 


45 


9 


64 


2 


57 


82 





21 






The table shows the properties for the reduced cluster with the number of protostars Afgink, the total cluster mass Mc, the average 
protostellar mass (M), the radius i?c, and the protostellar number density n^""^. Column ti indicates the time of the formation of the first 
protostar, Ati2 the time difference between the formation of the first and the second protostar. craD and ctid show the stellar velocity 
dispersion of the cluster. The key values for the mass segregation are the relaxation time trclaxi the available lifetime of the cluster (avail 
in units of the relaxation time, and the total and normalised number of protostars that had enough time to relax dynamically A^scg and 
fseg = .'Vseg/A^aink- Tho column "seg?" indicates the segregation state of the cluster: significantly mass segregated (H — h), marginally mass 
segregated (+), not mass segregated (0), and inversely mass segregated (— ). 
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Figure 9. Subclusters in the TH runs. The left picture shows TH-m-1 with the subclusters SCI and SC2. The two largest subclusters 
in TH-m-2 on the right are labelled SC3 and SC4. The circles indicating the subclusters' diameter are to scale. The total size of the plot 
is 0.13 pc in both x and y direction. 



4.3.2 BE runs 

The effects of the much more dominant central infall during 
the collapse of the BE setups can be seen in the average dis- 
tance between the sink particles and the Q- value in figure [TT| 



The separation distribution shows only one significant max- 
imum for all simulations. However, the peak for the BE-c-2 
run is at a much larger distance. There, the sink particles 
form along large elongated filaments and lead to larger mean 
separations than in the other BE setups. Here, the strong 
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Figure 10. Global cluster values for the TH runs. Left plot: TH-m-1 shows three different peaks in the distribution function (see equa- 
tion |ll| l: the one at 9,000 AU corresponds to the distance of SC2 to SCI (see figure[9]l, the peak at 13,000 AU to the degenerate distance 
of SC2 to SC5 and SC6 and the last peak describes the distance from SCI to SC5 and SC6, which is also degenerate within the width of 
the distance bin. TH-m-2 shows two main subclusters, whose distance corresponds to the peak at 15,000 AU. Right plot: After roughly 
100 sink particles have formed, the Q value (see equation |10| approaches a constant value which is similar for both of the runs, indicating 
a high degree of substructure in both clouds. 



effects of the compressive turbulent motions have a major 
impact. The mean separation for both runs with compres- 
sive turbulence is significantly larger than for the other runs 
(see (s) in table [3|. The Q values and the resulting degree 
of substructure are very different and strongly change with 
time (and consequently A'^sink) depending on where the sink 
particles form. BE-c-2 shows strong substructure from the 
very beginning, BE-c-1 forms protostars at larger radii at a 
later stage in the simulation, leading to a decrease of Q at 
around A^sink ~ 170. The two runs with the lower number of 
sink particles (BE-s-1 and BE-m-1) have the highest values, 
revealing a rather smooth cluster without much substruc- 
ture. 



4.3.3 PL15 runs 

The even stronger mass concentration in the PL 15 profiles 
shows a systematic influence on the mean distance between 
the sink particles. The mean particle separation for the 
PL15-C-1, PL15-m-2 and PL15-S-2 runs is roughly 15 - 35% 
smaller than in the corresponding BE runs (seetablejsf. The 
fact that the mean separation in PL15-C-2 is larger than in 
BE-c-2 is just due to the fact that the former one forms 
fewer sink particles; the positions of the distant sink parti- 
cles at large radii are similar. Figure[l2]shows the separation 
function and the Q values. The distribution function on the 
left shows one main peak for all setups. The peak for PL15- 
c-2 is much wider, reflecting the larger central cluster. In 
addition, the setup forms more protostars further out than 
other setups. In combination with the lower total number 
of particles than in the BE-c-2 case, this yields the large 
value of (s) and result in the lowest Q value for PL15-C-2. 
PL15-C-1 and PL15-m-2 are around the threshold value to 
substructure [Q = 0.8), PL15-S-2 is smooth over almost all 
the simulated time. 



4-3.4 Comparison 

There are some general trends of the subclustering proper- 
ties. The flatter the initial density profile is, the more im- 
pact has the turbulent velocity field. This causes collapsing 
regions to form at larger separations from each other. The 
observed relation (Qth) (Qbe) J; (Qplis) supports this 
intuitive picture. In a similar manner, compressive turbulent 
modes lead to collapsing filaments more quickly, not allowing 
the gas to assemble as close to the centre as in solenoidal tur- 
bulent cases. Therefore, within one density profile, the im- 
pact of turbulent modes shows (Qcomp) (Qmix) (Qsoi). 

4.4 Reduced Cluster Properties 

Having analysed the total set of protostars in the entire 
cloud, we now focus on the central regions of the main clus- 
ters in each setup, ignoring the outliers that do not belong 
to the main cluster. In order to find the individual compact 
clusters, we iteratively exclude outlier protostars until we 
reach a converged cluster configuration. We first select the 
main region by eye. In the two TH runs we select the already 
mentioned subclusters (see figure|9|, in all other setups with 
many sink particles we chose the central cluster. The par- 
ticle reduction method works as follows. We find the centre 
of mass of the set of particles. Then we compute the aver- 
age separation (s) between protostars and remove all objects 
that are located at radii larger than three times the mean 
separation from the centre of mass. We then recalculate the 
centre of mass and repeat the exclusion until no further par- 
ticle is excluded from the set of objects. The radius of the 
cluster Rc is set to 3 (s), ensuring that all selected particles 
are within the cluster radius. The factor three is somewhat 
arbitrary, but after some tests it turned out to be a useful 
distance factor that does exclude all very distant particles, 
but no or very few particles that could be dynamically im- 
portant for the cluster within the simulated time. The key 
values for the reduced clusters are listed in table |4l For the 
following discussion we focus on the reduced clusters. 

As the motions in the forming cluster are highly chaotic 
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Figure 11. Cluster properties for the protostars in the BE runs. The plot of the separations (left figure) clearly shows the formation of 
only one main cluster for all runs, indicated by only one main peak in the distribution of protostellar separations. However, the cluster 
structure varies significantly (right figure). The Q value differs by a factor of more than 5 for the individual runs and shows a correlation 
with the turbulent modes. Compressive modes show more substructure than mixed and solenoidal modes. 
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Figure 12. Cluster properties for the protostars in the PL15 runs. The plot of the separations (loft figure) clearly shows the formation of 
only one main cluster for all runs. However, the cluster structure varies significantly (right figure). 



and the number of protostars is constantly growing, the time 
evolution of the reduced cluster properties fluctuates, i.e., 
every time step, the reduction algorithm chooses different 
protostars to belong to the reduced cluster. It is therefore 
impossible to follow single protostars within the reduced 
clusters. In the further analysis we thus concentrate on the 
clusters at the end of the simulation. 



4.5 Mass Segregation 

We address the mass segregation problem in two ways. 
Firstly, we investigate the time that each protostar had for 
dynamical mass segregation after its formation, and sec- 
ondly, we analyse the reduced cluster at the end of the sim- 
ulation with the minimal spanning tree, neither taking into 
account the different formation times of the particles nor the 
change in mass during the accretion process. 

Although the degree of mass segregation can not be cal- 
culated for a single particle but has to be seen as a global 
cluster property, we analyse the possibility to dynamically 
mass segregate via two-body relaxation for every single pro- 
tostar. According to equation ( 16 1 we set the time tscg to 



the time that the sink particle had for mass segregation, 
i.e., the difference between the end of the simulation and 
the creation time of the protostar in question. From that 
we infer the threshold mass Alsag, with the given final val- 
ues of Rc and a. If the mass of this particular protostar is 
larger than the threshold mass, we count it for possible mass 
segregation. The quantity A^scg in table |4] refers to the total 
number of possibly mass-segregated objects; /scg denotes the 
fraction Nscg/Nc- The strong dynamical effects during the 
formation of the cluster result in significantly varying values 
for Rc, Nc, and a. However, the combined quantity in equa- 



tion ( 16 1 differs much less and serves as a remarkably stable 



estimate. With a roughly constant formation rate of proto- 
stars, a strong correlation between the protostellar number 
density (Nc/Rc) and /sog as found in the simulated clus- 



ters is not surprising (see figure 14 1. The segregation fraction 
/scg covers a very large range (0.03 — 0.85), indicating that in 
some setups almost all objects had enough time to dynam- 
ically mass segregate, while in others hardly any protostar 
can relax in the cluster. 

For the second approach, we analyse the mass segrega- 



tion at the end of the simulation according to equation ( 14 1 
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Figure 13. Velocity dispersion of the selected main (sub)clusters 
as a function of stellar density. The data points represent the 
clusters at the end of the simulation. The clusters show a weak 
correlation with a significant scatter. 
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Figure 14. Possible fraction of dynamically mass-segregated 
stars as a function of protostellar number density of the clus- 
ter. The data points represent the reduced clusters at the end 
of the simulation and show a strong correlation with only little 
scatter. 



The values for Amst as a function of A'^mst for all clusters 
are plotted in figure |15[ In order to keep the plots read- 
able, most of the curves are shown without errorbars; we 
included errorbars for the lowest curves that still differ from 
unity within a 1-a error in order to give some indication of 
the uncertainties involved. In order for mass segregation to 
be eminent, Amst needs to be significantly above unity for 
mass segregation and significantly below unity for inverse 
mass segregation. The upper panel shows Amst for the TH 
runs. All subclusters except for SCI show mass segregation 
up to at least A''mst ~ 30, i.e., the 30 most massive proto- 
stars form a compact subset of the cluster members around 
the centre of the cluster. Including more than the 30 most 
massive objects to the subset enlarges the spatial extent such 
that the position of the chosen subset is hardly distinguish- 
able from a random subset of the same number of cluster 
members. SC3 and SC4 show a significantly higher degree 
of mass segregation below A^mst ~ 20 and A'mst ~ 12, re- 
spectively. This corresponds to a minimum segregated sink 
mass of 0.074 M© in SC3 and 0.11 M© in SC4 and con- 



tains roughly 40% and 37% of the total cluster mass. Even 
higher values for Amst can be found in the BE setups (mid- 
dle panel). Here the central clusters in BE-m-1 and BE-s-1 
show Amst > 1.5 below A''mst ~ 35 — 45 and A^mst ~ 20, 
respectively. The minimum segregated mass in BE-m-1 is 



Me, 



0.1 Mq, the total confined mass down to this mass 



is about 75%, in the latter case Mscg = 0.17 Mq, containing 
around 40% of the cluster mass. If one includes the second 
bump of Amst between 20 < A''mst < 40 in BE-s-1, the 
measured contained mass that is segregated is roughly 58%. 
Among the PL15 density profile only one cluster shows sig- 
nificant mass segregation, PL15-C-1. Amst is greater than 
1.5 for TVmst 19- This gives a minimum segregated mass 
of Afseg = 0.11 Mq and corresponds to a fraction of about 
72% of the cluster mass. 

There is a weak correlation between the actual mea- 
sured mass segregation and the theoretically possible frac- 
tion of segregated protostars (/seg). The actual number of 
segregated stars iVMST.max with AMST(iV ^ A^MST.max) > 
1.5 is lower in almost all cases, but follows a consistent trend 
with A'sog. Taking into account that the protostars form at 
different positions and need some time to dynamically relax 
within the cluster, the relation A^MST,max < A^seg seems rea- 
sonable in comparison to an initially spherical cluster with 
a constant number of members. 

The actual mass segregation can also be compared to 
the total time that the cluster as a whole has for mass seg- 
regation. As the number of protostars changes with time, 
we count the available time starting at the point where two 
sink particles are formed until the end of the simulation. 
The ratio iavaiiAroiax in table |4] indicates how many mass 
segregation times the cluster evolved, again assuming that 
trciax at the end of the simulation is representative for the 
entire cluster evolution. There is again a weak correlation 
between this ratio and the degree of mass segregation. 



5 DISCUSSION 

In all simulations we set up cores with a very low ratio of 
kinetic to gravitational energy, i.e., the clouds are strongly 
bound. As the cores are isolated, they are disconnected from 
any potential dynamical impact from the surrounding en- 
vironment. The initially imposed supersonic turbulent mo- 
tions result in a global velocity dispersion for the gas of 
~ 0.5kms~^. Given the diameter of the core, 0.2 pc, this 
is close to the velocity dispersion we expect from Larson's 



relation (|Larson||1981 Solomon et al. 


1987 


Ossenkopf & 


Mac Low||2002| |Heyer & Brunt||2004| 


loman-Duval et al. 


2011 


I. However, it is lower than the observed turbulent ve- 



locity component of the massive dense cores in Cygnus X 
( |Bontemps et al.|[20To{ |Csengeri et al.||2011[ ). The observed 
cores with very similar key properties to our cores, i.e., 
mass, size, and temperature, show velocity dispersions from 
~ 0.5 — 3.5kms~^, higher than the turbulent velocity dis- 
persions in our numerical setups. Observations of massive. 



dense filaments show supersonic infall motions (Schneider 
|et al.|2010[ ), which may easily lead to a much more dynam- 
ical formation of the cores. We do not impose an initial net 
rotation to the core, however, the random turbulent pattern 
of high and low- velocity regions in different density environ- 
ments results in a net rotation of the cores with a ratio of 
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Figure 15. Mass segregation Amst fof a-U setups at the end of 
tlie simulation. For tlie setups wliere the deviation from unity is 
not obvious, Amst is plotted with errorbars. 



rotational to gravitational energy ranging from 10""'-10"^ 
in agreement with the values for the dense cores in Cygnus X 
( Bontemps et al.|2010 Csengeri et al. 2011[ ). During the sim- 
ulation, the velocity dispersion increases significantly due 
to the strong global infall and reaches values that are more 
consistent with the observed ones. After some ~ 10 — 40 kyr, 
depending on the initial density profile, the cores as a whole 
reach or exceed a virialised energy budget -Ekin/I-Epotl 5S 0.5. 
The final energy balance is in agreement with the theoreti- 



cal virial analysis in Shetty et al. (20101. They investigated 



the scaling relations between mass, size, and virial state of 
clumps of different sizes that formed self-consistently in tur- 
bulent flows. The virial state of their clumps with similar 
sizes and masses to our setups is consistent with our en- 
ergy analysis. Also the measured line widths of our cores 



is consistent with the analysis in Shetty et al. (20101. The 



et al. (20111 notice small-scale turbulent motions with high 
velocities (a few kms~^) in high-resolution studies of the 
central region of the cores. 

As soon as protostars form, the question of early sub- 
structure and mass segregation arises. These two properties 
of young stellar clusters can not be disentangled and anal- 
ysed separately. In particular, the determination whether a 
cluster shows primordial or dynamical mass segregation sen- 
sitively depends on the definition of mass segregation and 
spatial demarcation of the region in question. 

The theoretical analysis of a self-gravitating TV-body 
system predicts dynamical mass segregation via two-body 
relaxation and dynamical friction that an object experiences 
while moving through a sea of other objects. For different 
properties of the cluster, the dynamical friction and the re- 
sulting dynamical relaxation time of the total cluster differs 
(Chandrasekha r 1943; McMillan fc Portcgics Zwart 2003] 
Spinnato et al. 2003; Fc Uhauer fc Lin||2007[ ). The global re- 



laxation time, defined as a statistical quantity with only 
global cluster properties and thus not reflecting any sub- 
structure, therefore only serves as a rough estimate. Depend- 
ing on how well these global quantities fit the observed or 
simulated system, the relaxation time might differ signifi- 
cantly from the time scale of local dynamical interactions. 

The question whether dynamical mass segregation can 
be excluded based on a time-scale argument, can therefore 
only be answered for a specific definition of mass segregation 
and for a well-defined cluster or subcluster region. Tradition- 
ally, numerical work started without initial mass segregation 
and investigated the purely dynamical aspect of the A''-body 
system, without taking into account the dynamical changes 
of the individual A''-body objects, like mass accretion in the 
early phase of cluster formation and the mass loss due to 
winds. Recently, several p rescriptions of initia l mass segrcga- 
tion have been developed ( Baumgardt et al. 2008. , Subr et al.| 



|2008[ [Vesperini et al.|2009[ ), still investigating the cluster as 
a whole without local substructure. 

One basic problem with the analysis of mass segregation 
is the definition of what mass segregation actually means. 



Allison et al. ( 2009 1 use the minimal spanning tree (MST) of 



the most massive stars in comparison to the MST of random 
stars and thus define mass segregation as the most massive 
stars being located closer to each other than the same num- 
ber of randomly picked stars. As long as a single cluster 
or a conglomeration of several individual clusters does not 
show massive outliers, this methods works stably. In case of 
massive outliers, this method needs to be slightly modified 
( Maschberger fc Clarkel|2011[ [Olczak et al.||2011[ ). In obser- 
vational studies, mass segregation is mostly defined as more 
massive stars being located closer to the centre of the clus- 



ter (e.g. Hillenbrand 1997 Hillenbrand fc Hartmann||1998 



IFischer et al.'1998' 'de Grijs et al.'2002', 'Sirianni et al .]|2002 
Gouli ermis et al. 2004 , Hufl^ fc Stahlcr 2006, Stohe et al 
2006 r |Sabbi et al.i2008| |Gennaro et al.|2011| |Kirk fc Myers' 



2011 1. However, the definition of the centre of a young star 



increasing values for aio are dominated by the gas motions 
in the dense central region, which is also observed. ICsengerT] 



forming region with a large degree of substructure is not 
obvious. 

One possibility to study mass segregation in resolved 
clusters is to investigate radial variations of the initial stellar 
mass function. In unresolved clusters the different inferred 
radii in different wavelengths may indicate mass segregation. 
However, in both cases, mass segregation is difficult to iden- 
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tify given the observational difficulties (e.g., Ascenso et al, 
2009[ [Portegies Zwart et al.|2010[ |. 



Even more difficult is the answer to the question about 
primordial versus dynamical mass segregation. In order for 
global mass segregation to be primordial in nature, it is re- 
quired that stars with a given mass m must be more cen- 
trally concentrated than stars with the average stellar mass 
(m) and that the cluster must be younger than the dynam- 
ical friction time scale for that given mass m, i.e., the more 
massive stars must have formed closer to the centre. This 
global picture is consistent with numerical simulations (e.g., 
Klessen & Burkert||2000l iBonnell & Bate 2006). However, 



this time scale argument only holds for spherical clusters 
in virial equilibrium. If clusters form through mergers of 
smaller subclusters, these subclusters might have enough 
time to dynamically relax and mass segregate because of 
the much smaller size and the higher number of stellar en- 
counters. The degree of mass segregation in merged clusters 
is significantly higher than would be expected from a global 
time scale analysis ( McMillan et al.]|2007 Moeckel & Bon 
|nelll'2009). In addition lAllison et al I \20W\ show that dy 



namical mass segregation is very fast even without mergers 
of partially mass segregated substructures. Therefore, a de- 
tailed analysis of the formation of substructure and the col- 
lapse of stars within them is crucial to fully understand the 
mass segregation process. The analysis of our reduced clus- 
ters and subclusters with their own dynamical and orbital 
centre shows that there is a weak correlation between the 
possible degree of segregation /scg and the actual mass seg- 
regation. Given the fact that we follow the evolution of our 
clusters for only a very short time, it seems very likely that 
dynamical mass segregation can provide a significant contri- 
bution to the mass segregation within the subclusters. If, in 
addition, the bigger stellar cluster that formed by mergers of 
smaller subsystems, can inherit a reasonable degree of mass 
segregation of the progenitors, it becomes very difficult to 
exclude dynamical effects on different spatial and dynamical 
levels to be responsible for mass segregation of a cluster. 

A further complication in the mass segregation process 
arises from dynamical effects due to ejected stars. |Yu et alT] 
(20111 showed that removing ejected stars has an effect on 



measuring the mass segregation of the cluster. Likewise, the 
initial velocity distribution influences the segregation pro- 
cess. In our analysis, we do not take into account the ef- 
fects of escaping stars. In fact, we do not have stars that 
entirely escape from the cloud. Whether this is due to the 
low number of objects in the central region of the cluster, 
the relatively short evolutionary time of the simulation or 
due to the gravitational softening, remains an open ques- 
tion. Follow-up simulations that evolve the clusters for a 
longer time and with accurate protostellar sizes and the re- 
sulting gravitational potential are needed to clarify this ef- 
fect. Concerning the initial velocity distribution, our setups 
show significant differences from the simulations by |Yu et al.| 
( |2011[ ). As we follow the formation of protostars in the gas 
cloud, the initial protostellar velocity distribution is not a 
free parameter, but is inherited from the gas motions of col- 
lapsing regions. In addition, the protostars in our clusters 
are embedded in a dense cloud whereas the simulations by 



Yu et al. (20111 only consider the motions of the particles 



The total cluster including all protostars shows a sub- 
virial energy budget, indicating that the relaxation time is 
larger and thus the dynamical mass segregation process of 
the total cloud is slower than in a virialised case. How- 
ever, the central regions, where the crossing times are much 
smaller and stellar encounters more frequent, the iV-body 
system is virialised. The central region therefore does not 
suffer from a dynamical delay concerning the mass segrega- 
tion process. In addition, the simple analysis of dynamical 
mass segregation does not include the effects of gas, but only 
the dynamical friction due to the other stellar objects in the 
sample. In addition, the gas also provides dynamical fric- 



tion fDokuchacv"1964"Rudcrman & Spicgcri971 Rephaeli 



& Salpeter. ,1980^ Pstriker J999, ^Lee 



Stahler||2011[ ). Due 



to the turbulent motions, an analytic estimate is difficult to 
apply in our collapsing core. Nevertheless, this additional 
friction helps to increase the dynamical cross sections and 
thus makes stellar encounters more frequent, resulting in an 
acceleration of the dynamical mass segregation. 

An interesting aspect that weakens the effect of dynami- 
cal mass segregation is presented in recent work by |Converse| 
|fc Stabler (2011), where they argue that low-A^ systems with 
an even higher number of objects than in our clusters do 
not relax dynamically. If this also applies to accreting stel- 
lar systems with gaseous background, a large degree of mass 
segregation might not be possible in the smallest subclusters 
but only later after some merger events. We nevertheless do 
not expect dynamical relaxation to become completely irrel- 
evant because of the low number of protostars in our clusters 
and subclusters. 

As a remark, we want to point to recent studies by |Krui-] 
jssen et al. (20111. They analysed the substructure within 



clusters as well as the dynamical state of the stellar clus- 
ter when gas expulsion becomes important, i.e., at a slightly 
later stage of the evolution of the cluster. Analysing the 



simulations of Bonnell et al. ( 2003 2008 1, they find that the 



stellar system quickly reaches a globally virialised state if 
the gas potential is excluded and the stellar system is fol- 
lowed with pure A'^-body dynamics. Their results support 
the evolutionary picture of the formation of protostars that 
we see in our simulations. New protostars that form at larger 
radii from the centre of the cluster in gas dominated regions 
have sub- virial velocities. As soon as they decouple from the 
gas motion and move to the central gas-poor environment, 
they quickly virialise. The analysis of the simulations by 
Bonnell et al. ( 2003 2008 1 with a focus on mass segregation 



(Maschberger & Clarke 20111 shows global mass segrega- 



tion from very early times which continues throughout the 
simulation. This is also in agreement with our results. Fur- 
thermore, the degree of mass segregation is only mildly infiu- 
enced during subcluster merging, which suggests that early 
mass segregation can survive strong dynamical impacts. 

Note that the number density of protostars in the cen- 
tral region of the clusters is high enough for protostellar col- 
lisions to become important ( Baumgardt fc Klessen|["2011 1. 
This could indeed lead to changes in the stellar initial mass 
function. A discussion of the initial stellar mass function for 



all setups is presented in Girichidis et al. (2011a I 



without background gas. Finally, their simulated times are 
orders of magnitude longer than in our case. 



Concerning the physical processes, we chose a simple 
setup neglecting radiative feedback, magnetic fields, jets and 
outflows from the young protostars, and chemical reactions. 
Previous studies have shown that magnetic fields tend to 
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reduce the fragmentation (|Ziegler 


2005 Banerjee & Pudritz 


2006| 


Hennebelle & Teyssier 2008 


Wang et al.|2010 


Biirzle 


et al. 


20111 


Hennebelle et al. 201 1| Peters et al. 


2011a| 


Seifried et al. 


2011 1 without preventing it entirely. We there- 



fore expect our setups to form fewer stars in a magnetised en- 
vironment. Similar effects would be expected if we included 
radiative feedback. |Bate| poMl , IKrumholz et al.l \2009\ , |Pe" 



ters et al. 



( 2010a|( 



b I found reduced fragmentation in sim- 



ulations, without suppressing it entirely. |Commergon et aT] 
(20111 and Peters et al. (2011b I combined both, magnetic 



fields and radiative feedback, finding that the complex inter- 
play between the two processes reduces fragmentation with- 
out entirely suppressing it. Concerning mechanical feedback. 
Dale & Bonnell (2008) find that winds from massive stars 



can slow down the star formation process, but that the time- 
scale on which they can expel significant quantities of mass 
from the cluster is of the order of 10 free-fall times, which 
is much longer than the simulated time of our clusters. As 
their cores are not as dense and unstable as ours, we expect 
the effects of winds to be much less significant. A difi'erent 
fragmentation behaviour might well show differences in the 
formation pattern of stars, see also the different formation 



modes in Girichidis et al. (201 lb I. The substructure within 



the clusters is also expected to be influenced by different 
physical processes. In contrast, the overall energetics seems 
to be dominated by the global gravitational collapse of the 
cloud. 



6 SUMMARY AND CONCLUSION 



We analysed the simulations described in [Girichidis et al.] 
( 2011a|b l with the focus on the properties of the embedded 
young stellar clusters. We analysed the energy evolution of 
the gas and the nascent cluster, computed the degree of 
subclustering, and quantified the mass segregation in the 
continuously growing clusters. Our main conclusions can be 
summarised as follows. 

In all setups, the collapsing cloud virialises within the 
simulated time, which corresponds to a star formation ef- 
ficiency of 20%. Just considering the gas, all clouds have 
a virial or super- virial energy budget Skin > 0.5|-Epot|, the 
runs with only one protostar have significantly higher ra- 
tios of kinetic to gravitational energy. Although the total 
mass of all protostars is only 20% of the total cloud mass, 
their total kinetic energy is larger than that of the gas in 
the cases with multiple protostars. In contrast, the three 
runs with only one protostar show a smaller ratio of ki- 
netic energy of the protostar to kinetic energy of the gas, 
which can be explained by the vanishing momentum impact 
of opposite accretion fiows. Analysing the entire stellar clus- 
ters as pure A'^-body systems, we find an overall sub-virial 
energy balance with _Ekin ~ 0.2|_E'pot|, independent of the 
varied initial conditions. If we concentrate on the central 
regions of the clusters (innermost ~ 10 — 30% of the proto- 
stars), we find virialised conditions. This difference can be 
explained by the formation history of the cluster. New pro- 
tostars continue forming at increasing radii from the centre 
of the cloud due to the lack of available gas in the central 
region ( Girichidis et al.|2011b K These protostars inherit the 
kinetic energy from their parental gas region, which is rela- 
tively low in comparison to their gravitational contribution. 



i.e., new stars form at sub-virial velocities. Soon after their 
formation, the protostars decouple from the gas and agglom- 
erate in the central region, where they virialise. 

The degree of subclustering strongly depends on the ini- 
tial density profile. Initially uniform density allows for tur- 
bulent motions to form distinct subclusters before the global 
collapse can confine the gas in one cluster in the central re- 
gion. With a Q value of ~ 0.2, these clouds show consid- 
erable substructure with distinct conglomeration of proto- 
stars. The stronger the initial mass concentration around 
the centre of the cloud, the less subclustering is found. 
Bonnor-Ebert-like spheres show mainly one dominant cen- 
tral cluster with some substructure. The considered power- 
law density distributions form more compact protostellar 
clusters with less internal structure, if they form clusters 
at all. In three strongly condensed setups the cloud does 
not fragment and forms only one protostar. In general, we 
find that the Q parameter, used to quantify subcluster- 
ing, shows the following trend: (Qth) < (Qbe) < (Qplis), 
where lower Q means more substructure. We also note dif- 
ferent subclustering trends with different turbulent modes. 
For a given density profile, compressive modes lead to a 
higher degree of substructure than mixed modes, which in 
turn lead to more substructure than solenoidal modes, i.e., 

(Qcomp) (Qmix) ^ (Qsol). 

Focusing on the central region of the clusters, where 
outliers are removed from the set of protostars, roughly half 
of the clusters show mass segregation. The degree of mass 
segregation varies strongly between the clusters, however, 
no cluster with significant inverse mass segregation is found. 
Except for one cluster (PL15-m-2), the mass segregation ra- 
tio does not drop below 0.5. The mass segregation is con- 
sistent with the time for dynamical mass segregation, so all 
the clusters had enough time for dynamical relaxation of the 
most massive objects in the cluster. In the simulated col- 
lapsing cores, primordial mass segregation is not necessarily 
required to achieve a significant mass segregation at the end 
of the simulation. However, due to the ongoing formation of 
protostars and the increase in protostellar mass due to ac- 
cretion, the cluster is exposed to continuous momentum and 
energy impact from the surrounding gas, which may mod- 
ify the actual mass segregation behaviour in comparison to 
the idealised process of dynamical mass segregation via two- 
body relaxation. A contribution that may have a significant 
influence is the episodic accretion of gas as well as the fact 
that the protostars follow the global flow pattern of the gas 
they form from, before they dynamically decouple from the 
gas. Overall, there is no clear correlation between the ini- 
tial conditions and the mass segregation in our simulated 
clusters. 

We conclude that the kinetics of young stellar clus- 
ters do not strongly depend on the initial density profile, 
nor on the initial structure of the turbulent modes. This 
is because the nascent protostars quickly decouple dynami- 
cally from the parental filament in which they were formed. 
The interactions as an A'^-body system dominate the clus- 
ter motions. Continuous formation of subsequent protostars 
with initially sub-virial velocities lead to a globally sub-virial 
(i?kin/|-Bpotl < 0.5) state for the majority of the protostars. 
Taken into account the dynamics of small subclusters with 
dynamical times much smaller than the dynamical time of 
the entire cloud, the measured degree of mass segregation 
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is fully consistent with dynamical mass segregation, there is 
no need for primordial mass segregation in our simulations. 



APPENDIX A: GRAVITATIONAL FORCE 
SOFTENING 

We used the gravitational softening for the sink particles as 
described in Price & Monaghan ( 2007 1. The potential energy 

(Al) 



can be written as 



-Bpot = '^Gm.imj (f>{ri - rj,h), 

where h is the smoothing length, which is set to the accretion 
radius of the sink particles h = raccr/2, and (j)(r, h) is given 
by {q = r/h) 



4>{r,h) 



h- 



1 /• 2 „2 3 „4 I 1 „5 7\ 



i?^ - is?" + To^^ - i) . 



^ g < 1 

1 s: g < 2 

(A2) 

Note that the function (j> is defined such that the potential 
energy is negative. 
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